Hidden symmetries of two-electron quantum dots in a magnetic field 



N. S. Simonovic 1 and R.G. Nazmitdinov 1,2,3 
1 Max- Planck- Institut fur Physik komplexer Systeme, D-01187 Dresden, Germany 
2 Departament de Fisica, Universitat de les Illes Balears, E-07071 Palma de Mallorca, Spain 
'tubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia 



Using a classical and quantum mechanical analysis, we 
show that the magnetic field gives rise to dynamical symme- 
tries of a three-dimensional axially symmetric two-electron 
quantum dot with a parabolic confinement. These symme- 
tries manifest themselves as near-degeneracies in the quan- 
tum spectrum at specific values of the magnetic field and are 
robust at any strength of the electron-electron interaction. 
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A three-dimensional harmonic oscillator with frequen- 
cies in rational ratios (RHO) and a Coulomb system are 
benchmarks for the hidden symmetries which account for 
the accidental degeneracies of their quantum spectra (see, 
e.g., [1]). Recent advances in nanotechnology create a re- 
markable opportunity to trace their dynamical interplay 
in mesoscopic systems like quantum dots, which confine 
a few electrons to a space of a few hundred nanome- 
ters. Indeed, competition between a confining poten- 
tial, approximated quite well by the harmonic oscilla- 
tor (HO), and the repulsive electron-electron interaction 
produces a rich variety of phenomena, for example, in a 
two-electron quantum dot (QD) under a perpendicular 
magnetic field (see [3,2] and references therein). In fact, 
a two-electron QD becomes a testing-ground for differ- 
ent quantum-mechanical approaches [2] and experimen- 
tal techniques that could provide highly accurate data 
for this system [4,3]. 

If the HO and the Coulomb potential are combined, 
most of the symmetries are expected to be broken. Nev- 
ertheless, in particular cases, the Coulomb (Kepler) sys- 



tem and the RHO may have common symmetries, as it 
was already noticed a long time ago [5]. The authors of 
Ref. [5] could not find, however, a physical application 
for this phenomenon. These symmetries were rediscov- 
ered in the analysis of laser-cooled ions in a Paul trap 
[6] and of the hydrogen atom in the generalized van der 
Waals potential [7]. The major aspect of the present pa- 
per is to demonstrate that the hidden symmetries could 
be observed in a two-electron QD under a tunable per- 
pendicular magnetic field if the effective confinement po- 
tential is indeed the three-dimensional (3D) HO. To this 
aim we focus our analysis upon the nonlinear classical 
dynamics of the system. At certain conditions the mo- 
tion becomes integrable and this indicates the existence 
of the symmetries in the quantum spectrum. 

The Hamiltonian of a two-electron QD in a magnetic 
field reads 

2 



H = J2 



ri - r 2 



(1) 



where a = e 2 j \-ksqE t . Here, e, m* , Eg and e r are the 
unit charge, effective electron mass, vacuum and rela- 
tive dielectric constants of a semiconductor, respectively. 
For the perpendicular magnetic field we choose the vec- 
tor potential with a gauge A = x r = ^B(—y, x, 0). 
The confining potential is approximated by a 3D axially- 
symmetric HO U(r) = m*[wQ (x 2 + y 2 ) +uj 2 z 2 ]/2, where 
Huj z and hujQ are the energy scales of confinement in the 
z-direction and in the xy-plane, respectively. In contrast 
to a 2D description of the QD this approximation pro- 
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vides a consistent description [8] of various experimental 
features [9] : the energy spectrum for small magnetic field, 
the value of the magnetic field for the first singlet-triplet 
transition, and the ratio of the lateral to vertical exten- 
sion of the dot. In the present analysis we neglect the 
spin interaction, since the corresponding energy is small 
compared to the confinement and the Coulomb energies. 

Introducing the relative and center-of-mass (CM) co- 
ordinates r = ri — r 2 , R = |(ri + r 2 ), the Hamilto- 
nian (I) is separated into the CM and relative-motion 
terms H = Hcm + Hrei [10]. The solution to the first 
term is well known [11]. It possesses all symmetries of 
the HO, since the Coulomb interaction enters the second 
term for the parabolic confinement [12]. In the follow- 
ing we concentrate on the dynamics of i? r oi- For our 
analysis it is convenient to use cylindrical scaled coordi- 
nates, p = p/l , p P = Pph/Ti, z = z/l , Pz = Pzh/K, 
where Iq — (h/ p.uja) 1 / 2 is the characteristic length of the 
confinement potential with the reduced mass p, — m*/2. 
The strength parameter a of the Coulomb repulsion goes 
over to A = 2a / {TiujqIq) . Although our consideration 
is general, for the demonstration we choose the values 
hujo w 2.8 meV and u z = 2.5 uj which are close to those 
obtained in our 3D analysis [8] of the experiment [9] . For 
the effective mass m* — 0.067m e and the dielectric con- 
stant s = 12, which are typical for GaAs, the value A = 3. 
Hereafter, for the sake of simplicity, we drop the tilde, i.e. 
for the scaled variables we use the same symbols as before 
scaling. 

In these variables the Hamiltonian for the relative mo- 
tion takes the form (in units of hcoo) 

, ^rcl 1 / 2 , m2 , 2, -22, -22, ^\ 

h= n^ = 2\ p p + -^ +p ' + u p p + U) * Z + r)- UJLm 

(2) 



where r = (p 2 + Z 2 ) 1 / 2 , iD = w/ljo, tti = l z /h. Here, ujl — 
eB/2m*c is the Larmor frequency and uj p = (lo 2 l + w 2 ) 1 / 2 
is the effective confinement frequency in the p-coordinatc, 
which depends through uul on the magnetic field. In this 
way the magnetic field can be used to control the effective 
lateral confinement frequency of the QD for a fixed value 
of the vertical confinement, i.e. the ratio lu z /lu p . 

Beside the energy (e = h) , the z-component of angular 
momentum l z is also an integral of motion due to the 
axial symmetry of the system. Therefore, the magnetic 
quantum number m is always a good quantum number. 
Since the Hamiltonian (2) is invariant under the reflection 
of the origin, the parity n is a good quantum number too. 

Although the motion in ip is separated from the mo- 
tion in the (p, z)-plane, the problem is in general non- 
integrable, since the Coulomb term couples the p and 
^-coordinates. Examination of the Poincare sections by 
varying the parameter uj z /uj p (see Fig. 1 for examples) in 
the interval (1/10,10) with a small step indicates that 
there are five integrable cases. The trivial cases are 
uj z /ujp — » and uj z /uj p — ► oo, which correspond to ID 
circular and 2D vertical QDs, respectively. The nontriv- 
ial cases are uj z /llJ p = 1/2, 1,2. These results hold for any 
strength of the Coulomb interaction and agree with the 
results for the Paul trap [6]. Below we discuss the non- 
trivial cases only. The typical trajectories in cylindrical 
coordinates are shown in Figs. 2a,c. 

The results obtained with the aid of the Poincare sur- 
faces of sections are invariant under the coordinate trans- 
formation. On the other hand, the integrability is a nec- 
essary condition for the existence of a coordinate system 
in which the motion can be separated. In turn, the analo- 
gous quantum mechanical system would be characterized 
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by a complete set of quantum numbers. 

At the value u' L = {ojI—lj 2 ) 1 / 2 the magnetic field gives 
rise to a spherical symmetry {lo z /lo p = 1) in the axially 
symmetric QD [8]. In this case the Hamiltonian (2) is 
separable in (scaled) spherical coordinates (see Eq. (7) 
in [8]). The additional integral of motion is the square 
of the total angular momentum l 2 . The spherical coor- 
dinates are a particular limit of the spheroidal (elliptic) 
coordinates well suitable for the analysis of the Coulomb 
systems (sec, e.g., [13]). Therefore, to search the separa- 
bility for the other integrable cases it is convenient to use 
the spheroidal coordinates (£,r],ip), where £ = (r 1 +r 2 )/d 
and r}= {r\—r 2 )/d. In the prolate spheroidal coordinates 
r i = [P 2 + (z + d) 2 ] 1 / 2 , r 2 = r. The parameter d € (0, oo) 
is the distance between two foci of the coordinate system 
(with the origin at one of them). In the limit d — > the 
motion is separated when uj z /ui p — 1 (Fig. 2b). In this 
limit £ — > oo such that r = d£/2 is finite, 77 = cos$, and 
we obtain the spherical coordinate system. 

Let us turn to the case uj z /uj p = 2 which occurs at the 
value of the magnetic field oj'[ = (w 2 /4 — w 2 ) 1 / 2 . In the 
prolate spheroidal coordinates the motion is separated in 
the limit d — > 00 (Fig. 2d). In fact, at d — > 00: £ — > 1, 
77 — > 1 such that £1 = d(£ — 1), £2 = d(l — r]) arc finite, 
- we obtain the parabolic coordinate system (£1 , £2 , f) 
where £1,2 — r±z. In these coordinates the Hamiltonian 
(2) has the form 



c = a z - uj 2 p 2 z 



(4) 



a + 6 
# (£ + *!) + A 



- ^ '—{h'h 



(3) 



and the equation (£1 + £ 2 )(h — e) = is separated into 
two decoupled equations for £1 and £2 variables. Simple 
manipulations define the separation constant 



which is a desired third integral of motion. Here a z 
is the z-componcnt of the Runge-Lenz vector a = p x 
1 + Ar/(2r), which is a constant of motion for the pure 
Coulomb system (i.e. when ui p = uj z = 0) [1]. The quan- 
tum mechanical counterpart of the integral of motion, 
Eq. (4), does not commute with the parity operator and 
we should expect the degeneracy of quantum levels. 

Due to the separability of the motion in the parabolic 
coordinate system, the eigenfunctions of the correspond- 
ing Schrodingcr equation can be expressed in the form 
ip(r) = /2(a) e lmip , where the functions fj are so- 

lutions of the equations 



dfj 



1 



+A - (-l)*2c 



fj=0, 3 = 1,2. (5) 



Eqs. (5) can be solved numerically and the eigenenergics 
and eigenvalues of c are determined iteratively by varying 
simultaneously e and c until the functions fj fulfill the 
boundary conditions: fj ~ £j m '^ 2 , ^jdfj/d^j = \m\fj/2 
for £j — > and fj — > for £j — > 00. Let n\ and n 2 be 
the nodal quantum numbers of the functions /1 and f 2 , 
respectively. Note that Eqs. (5) are coupled by the con- 
stants of motion and, therefore, both functions depend 
on all three quantum numbers (ni,n 2 , m). However, the 
simple product of these functions have no a definite par- 
ity. Since r -> -r {£1 -> £ 2 , 6 -> 6 > <fi -> <fi + the 
even/odd eigenfunctions are constructed as 



±( _ 1)m/ („ 1 , m)(a)/ („ 2 , m)(6)] 



(6) 



where N = rii+n 2 and k = \n\— n 2 \. These states are the 
eigenfunctions of h, l z , \c\ and the parity operator. For 
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|c| > the eigenstates (6) appear in doublets of different 
parity and, therefore, of a different total spin. For c = 
in Eqs. (5) /i = /2 and, obviously, only the states with 
parity 7r = (— l) m exist. 

For the magnetic field uj'[' = (Auj 2 z — uJq) 1 ^ 2 we obtain 
the ratio <jJ z /lu p = 1/2. The Hamiltonian (2) expressed in 
the oblate spheroidal coordinates (ri = [z 2 + (p + d) 2 } 1 / 2 , 
T2 = r) is separated for m = (at d — > oo). For m =/= 
the term m 2 / p 2 and, consequently, the Hamiltonian (2) 
is not separated in these coordinates. Note, for m = the 
cases uj z /ujp = 1/2 and 2 are equivalent if we exchange 
the p and z coordinates and, hence, the additional inte- 
gral of motion is \a p — Cj 2 z 2 p\. For m ^ we use the 
procedure described in Ref. [6] and obtain the following 
integral of motion 

C = [(a p u 2 z z 2 p) 2 + a% + 4m 2 ^r 2 ]V2, (7) 

where a p and a v are the p and ip components of the 
Runge-Lenz vector, respectively. Due to existence of 
three independent integrals of motion, h, m and c, which 
arc in involution, the dynamics for m/0, although non- 
separable, is integrable. The further analysis for m = 
is similar to the previous one and we omit it here. 

Let us denote by h* the Hamiltonian (2) for a specific 
value of the magnetic field when the system becomes sep- 
arable, i.e. for w£ = u' L , uo'[ or uj'[' (for m = 0). Then for 
an arbitrary value of ojl we can write 

h = h* + ^Au>ip 2 - Aw L m, (8) 

where A(i| = Co 2 — Cj* 2 , ACjl — & l — &l an d the term 
\ ACblp 2 is the only non-diagonal part of h in the eigenba- 
sis of h*. The eigenenergies of the Hamiltonian (2) (see 
Fig. 3) have been calculated with the use of the basis 
(6) and the spherical basis in the intervals < o>l < 1.5 



and the 1.5 < ujl < 5, respectively. The radial parts of 
the spherical eigenfunctions and f ni ,n 2 ,m, as well as the 
corresponding factors in the matrix elements (tpi \ p 2 \tpj), 
are evaluated numerically. The complete spectrum of the 
two-electron QD (Fig. 3a) shows the accumulation of lev- 
els with different quantum numbers into well pronounced 
bands at strong magnetic field. There is no obvious man- 
ifestation of the symmetries discussed above. In fact, 
the effects of symmetries are shown up for separated m- 
manifolds only (Fig. 3b). 

For non-interacting electrons (A = 0) , the energy levels 
of the QD arc Fock-Darwin levels [11] 

e = u) p (2n p + \m\ + 1) + Cj z (n z + \) - u) L m. (9) 

For rational ratios of uoz/uop the energy levels Eq. (9) 
are degenerate. It is simply the spectrum of the RHO 
in the external field Cj^m. For instance, at w z /w p = 2 
we have e = uj z (N + \m\/2 + 1) — uj'^m. The quantum 
number TV = n p + n z = ri\ + ni = 0, 1, 2, ... and each 
m-manifold consists of the shells characterized by this 
quantum number. Since the eigenenergies of the term 
ffcM with the corresponding quantum numbers are de- 
termined by Eq. (9) as well, the shells in the total spec- 
trum of the QD are not affected. 

The Coulomb interaction destroys the general symme- 
try of the 3D HO. However, the magnetic field can re- 
cover symmetries which are common for the RHO and 
Coulomb systems. At a relatively low value of the mag- 
netic field (for our parameters B k, 2.4 T) we re- 
veal the first manifestation of the hidden symmetries. 
This symmetry is determined by the integral of motion, 
Eq. (4) . It results in the appearance of shells at each m- 
manifold (Fig. 3b). There are exact crossings and repul- 
sions between levels of different and of the same parity, 



respectively, in each shell. The near-degeneracy of the 
quantum spectrum is reminiscent of a striking degener- 
acy observed for the RHO or pure Coulomb systems. At 
higher values of the magnetic field u' L (B w 7.5 T), the 
dynamical spherical symmetry appears, since l 2 becomes 
an additional integral of motion. This symmetry mani- 
fests itself as the attraction between levels with different 
orbital quantum numbers and the same parity (Fig. 3b). 
In contrast to spectra of pure Coulomb systems or of 
the RHO, there are no crossings between eigenstates of 
the subset characterized by a given quantum number m, 
since the accidental degeneracy is removed. Although 
the symmetry is recovered at very strong magnetic field 
(B ss 15.9 T) due to the appearance of the integral of 
motion Eq. (7), the dynamics is non-separable for m ^ 0. 
Note that shells are similar to the spherical case. 

The symmetries may be detected by studying the con- 
ductance of two-electron QDs at low temperatures. In 
particular, at uu'[ in the excited states there is the on- 
set of a singlet-triplet degeneracy related to crossings of 
the eigenstates (6) with |c| > (see Fig. 3b). The total 
spin S alternates between 1 and and the addition of a 
second electron with a spin-up or spin-down orientation 
to the QD will cost the same energy. At zero magnetic 
field B — two electrons occupy the same state with 
S = 0. At uj'I it becomes favorable for one electron to 
flip its spin. The electron reconfigures the charge and 
polarizes the two-electron QD leading to a Kondo type 
effect [14], i.e., to the increase of the conductance at low 
temperatures. The enhancement of the conductance in 
QDs due to the singlet-triplet degeneracy induced by the 
magnetic field have been discussed in [15-17] (for review 
see [18]). According to Rcf. [18], when the system is 



tuned to the degeneracy point the differential conduc- 
tance, dl/dV, would exhibit a peak at zero bias. The 
increase of the strength of the magnetic field removes 
the degeneracy and the peak will split onto two peaks re- 
flecting the single-particle spacing between singlet-triplet 
states. It should be noted, however, that this prediction 
is obtained in Ref. [18] within a schematic model where 
single-particle levels and the magnetic field are adjustable 
parameters. In addition, the electron-electron interac- 
tion is assumed to be weak. According to our analysis 
the onset of the singlet-triplet degeneracy holds at any 
strength of the electron-electron interaction in QDs. We 
suggest the mechanism, related to the hidden symme- 
tries, responsible for the occurrence of this degeneracy in 
two-electron 3D QDs at a certain value of the magnetic 
field. 

Summarizing, we have shown that the axially symmet- 
ric 3D quantum dot with two electrons exhibits hidden 
symmetries at certain values of the magnetic field. In 
particular, due to these symmetries the onset of a singlet- 
triplet degeneracy in excited states is found when the 
magnetic field value is uj'[. Finally, we are thankful to 
Jan-M. Rost for valuable discussions and constructive 
suggestions. 
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Figure Captions 



FIG. 1. Poincare surfaces of sections z — 0, p z > of the 
relative motion (A = 3, e = 10, m = 0) with: (a) lu z /lu p — 5/2, 
(b) LOz/ujp — 2 and (c) lo z /lu p = 3/2. The section (b) indicates 
that for the ratio lo z /lu p — 2 the system is integrable. 

FIG. 2. Typical trajectories (e = 10, m = 1) of the relative 
motion at A = 3 for lj z /lu p = 1 (a,b) and lo z /lu p — 2 (c,d) 
are shown in cylindrical and prolate spheroidal coordinates, 
respectively. 

FIG. 3. The lowest eigenenergies of Hamiltonian (2) (in 
units TwJo) at lo z /lo = 2.5 and A = 3 as functions of the ratio 
ojl/oJo for: (a) all m; (b) m = 0. The upper energy limit 
is chosen high enough to amplify the shell structure of the 
spectrum. The integrable cases: lu z /lu p — 2, 1 and 1/2 are 
indicated by vertical dotted lines (uil/wo = 0.75, V21/2 and 
\/24, respectively). 
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